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Abstract: We continue the study of the distribution of nonsupersymmetric flux vacua in 
lib string theory compactified on Calabi-Yau manifolds, as in hep-th/0404116. 
We show that the basic structure of this problem is that of finding eigenvectors of the 
matrix of second derivatives of the superpotential, and that many features of the results 
are determined by features of the generic ensemble of such matrices, the CI ensemble of 
Altland and Zirnbauer originating in mesoscopic physics. 

We study some simple examples in detail, exhibiting various factors which can favor low 
or high scale supersymmetry breaking. 
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1. Introduction 



In this work we continue the study of the statistics of flux vacua initiated in |18, ^, 19 1, 
and continue the study of supersymmetry breaking vacua begun in [^]. For a general 
introduction to this problem and its applications, and further references, see [^]. Some 
more recent works on vacuum statistics and related topics include |10, 24, |^. 

Our main result will be explicit formulae for the density of supersymmetry breaking 
vacua in an ensemble of effective supergravity field theories (EFT's) in which the superpo- 
tential is a "random variable" in a sense we define below. This type of ensemble includes 
the set of EFT's obtained by compactifying string/M theory with fluxes, in the limit that 
flux quantization can be neglected. We will give specific results for type lib flux compact- 
ification on Calabi-Yau |26], which at present is the best studied (and most calculable) 
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example, but stress that the techniques (and possibly many features of the results) are 
more general. 

Although we foresee many applications of such results to string duality and phe- 
nomenology, a question of primary current interest is the likelihood to discover super- 
symmetry at the energy range ITeV to be probed at LHC, perhaps as the traditional 
MSSM scenarios, "split super symmetry" ||2|, or perhaps other scenarios. As discussed 



in ||T^, 23], a systematic way to study this question is to try to count or estimate the 
number of string vacua which realize the various scenarios. Even with our limited present 
understanding, it is conceivable that we can argue that some scenarios are so unlikely that 
we should regard their observation as evidence against string theory, in other words we 
would have derived a strong prediction from string theory. Admittedly, this possibility 
depends on optimistic assumptions about the total number of string theory vacua, but in 
any case we believe vacuum counting will give essential information for any approach to 
string phenomenology. 

The main reason to reconsider traditional attitudes towards naturalness and super- 
symmetry is that, as becomes clear from a quantitative approach, the advantage gained 
by low energy supersymmetry is just not that large, compared to the current estimates of 
the number of string theory vacua; N^ac » 10^'^'^ vacua, based both on estimates of the 
number of flux vacua js], Q and on the idea that this large number of vacua is actually 
needed to solve the cosmological constant problem |^, 36, ^. 

To quantify the study of supersymmetry breaking, we adopt the following simplified 
description of the problem. We need to find the joint distribution of the Higgs mass, 
supersymmetry breaking scale and cosmological constant, 

among "otherwise phenomenologically acceptable" vacua, where 

Mlsy = Y.\PA? + Y.Dl. (1.1) 

A a. 

We then evaluate this at the (presumed) values Mh ~ lOOGeV ~ ^^"^"^ Mpianck and A ~ 
W~^^'^ Mpi^^^j^, and look at the resulting distribution of supersymmetry breaking scales. 
Rough arguments [^, |3^, ^] summarized in |^ suggest that for > MgusyMpianck, 
this goes as 

dl^lMn, Msusy, M ~ ,.4 "^^12 ir^^MMsusy]- (1-2) 

^"susyZ-^^planck ^^^planck 

To get a sense of what this means, consider a simple power law ansatz, 

d^l[Msusy\ ~ M^^sydMlsy 

Given such a distribution, the distribution of physical vacua would be weighed towards high 
scales if a > 2, showing that it would not take a huge growth in the number of vacua with 



breaking scale to dominate the advantage in solving the hierarchy problem. [35, 2^ Another 
way to make this point is to note that the ratio {Msusy Mgh/Msusy low)^ for two plausible 
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guesses at the supersymmetry breaking scale, the high scale Mgusy high ~ lO^^GeV and the 
intermediate scale Mgusy low ~ 10^'^GeV, is 10^^. This is not a large number in the present 
context. 

Of course the true distribution of supersymmetry breaking scales in string/M the- 
ory vacua will be far more complicated than this simple power law ansatz, with various 
components reflecting various supersymmetry breaking mechanisms and other structure 
in the problem. While it will clearly take a great deal of work to form any real under- 
standing of this true distribution, we will see in the following that certain generic features 
of the problem, following just from the the structure of supergravity and of generic dis- 
tributions of EFT's, do translate into specific properties and factors in the distributions, 
which might well be shared by the true distribution. Furthermore, many aspects of string 
compactification which at first sight look highly significant, such as mechanisms to pro- 
duce hierarchically small scales, can turn out to be no more important than these generic 
properties in the final statistics. 

We turn to detailed considerations of the problem of counting pure F breaking vacua 
in sections 2 and 3. Starting from a precise ensemble of EFT's such as those obtained by 
flux compactification, this is a mathematically well defined problem, but not a simple one. 
While we will obtain precise formulae as in 11 1, to some extent a more useful description 
of the basic result is that the distribution of supersymmetry breaking vacua is similar to 
that of super symmetric vacua, with certain "correction factors" which we will explain in 
detail. We also discuss the likely effects of taking flux quantization into account. 

In section 4 we make some comments on the distribution obtained by supersymmetry 
breaking by antibranes, or more generally by D terms. 

Section 5 contains conclusions. 



2. Preliminary considerations 

Our starting point is the M = 1 supergravity potential and its derivatives,^ 

V = e'^ {g°-^DaWDiW - 2,\wf^ + D'^ (2.1) 

daV = e'^{DaDbWD^W - 2DaWW) (2.2) 

DadbV = e'^iDaDbD.WD'W - DaD^WW) (2.3) 
D-ad^y = e'^iR^.abDdWD'^W + gMDcWD^W - D.WD-aW 

-2gf,-aWW + DbD,WD-aD'W). (2.4) 

For simplicity we will first assume the D terms are zero or at least independent of the fields 
(until section 4). 

^The covariant derivatives Da are both Kahler and metric covariant, i.e. when acting on W and its 
derivatives include the Kahler connection dalC, and when acting on tensors include the Levi-Civita connec- 
tion r^c- R is the curvature of the cotangent bundle, i.e. R^^^iXd = [Va, V^JXc — di{g'^'^dagce) Xd. At 
a vacuum dV — 0, it does not matter whether the outer derivatives are covariant derivatives or ordinary 
partial derivatives, because DdV = d^V. We set Mp = 1. 
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The approach we will take to finding statistics of vacua is to summarize the sum over 
all choices which go into the potential Eq. ( |2.1[) , in terms of a joint distribution dfi of the 
superpotential W and its derivatives, evaluated at a point z in the space of chiral field vevs, 
a candidate vacuum. To evaluate all of Eqs. ( |2.2[ )-(|2^) at a point z, the joint distribution 
must describe W{z), JC{z, z), and up to three derivatives of W at z. 

The explicit dependence of these quantities on /C(z), the Kahler potential at z, can be 
removed by either making redefinitions such as 1^ = e^^'^W (as in |11]), or equivalently 



making a Kahler- Weyl transformation to set IC{z,z) = 0. This does not change any of the 
equations because is covariantly constant with respect to Da- Let us now assume this 
has been done, and use the notations 

Fa = DaW{z); Zab = DaDbW{z); Uabc = DaDbDcW{z) (2.5) 

where capital indices A, B etc. are defined to be orthonormal complex indices with respect 
to the Kahler metric at z. We retain the symbol W for W{z). The tensors Z and U are 
symmetric in all indices. 

We can then rewrite Eqs. (|2.2D~(2.4) as 



SaV = ZabF^ - 2FaW (2.6) 

DaObV = UabcF"" - ZabW (2.7) 

DaObV = RabcdF^F^ + 8ab\F\'' - FaFb - 25ab\W\'' + ZacZbd- (2.8) 

Given a joint distribution 

d^i[W,FA,ZAB,UABci (2.9) 

such as might come from summing over all choices of flux in a given compactification, we 
will then compute densities such as 



p{z) = J dfi 6{V'{z))\detV"{z)\0{V") 
in terms of the joint distribution. 



2.1 Distributions 

The joint distribution Eq. ( |2.9D for an integral over fluxes in lib on CY in the large volume 
limit was worked out in To summarize, in these effective theories the Kahler potential 
/C is independent of the flux, while the GTVW formula 

W = J GAn{z) (2.10) 

tells us that the superpotential is linear in the flux G. Thus, W{z) and all of its derivatives 
are linear in the flux, and one can change variables in the integral over fluxes to find a 
simple joint distribution in which W, Fq, Fj and Zqj are independent. Here the index '0' 
refers to the dilaton, and / = 1, ... ,n to the complex structure moduli. Using threefold 
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special geometry, for the particular case of lib flux vacua, the variables Zab and Uabc are 
determined in terms of Z/ = Zq/ and Fa] for example 



Zij = TijkZ^ ■ 

Thus, the resulting distribution is 

= SW (f^+^FA (f'Zoj 6{L - \W\^ + - \Zoi\^) (2.11) 

with the other variables determined in terms of these. One could also express the distribu- 
tion for the original problem with quantized fluxes in this way, as a sum over lattice points 
embedded in the parameter space in a way determined by z and the periods. 

The superpotential Eq. ( |2.10D does not include Kahler moduli and as in |11| we will 



ignore these moduli. Of course, in an exact description, these must be taken into account. 
In particular, their stabilization requires a sufficient number of nonperturbative contri- 



butions to the superpotential |33, 12, 13], which may or may not be present in a given 
model. Moreover, to stabilize the volume at a reasonably large value, moderately small 
values of VFfiux ai'e needed. However, if thus stabilized, the Kahler sector has relatively 
little influence on the statistics and properties of vacua, because its contribution to the 
potential is exponentially small, and because by far the main degeneracy of vacua comes 
from the different choices of flux. In particular, inclusion of Kahler moduli typically only 
produces small shifts in the vacua found by ignoring them. This was made more precise in 
section 4.1, and jl^. The constraint of small Wfluxi needed to stabilize the model in a 
controlled (large radius) regime, merely reduces the number of suitable vacua by a factor 
l^fluxP and anyway this condition will be met automatically for the phenomenolog- 
ically most relevant vacua we will study in the following, namely those with A ~ and 

The requirement of a sufficient number of nonperturbative contributions is more subtle, 
but the results of |l^ suggest that many models satisfying this should exist, although 
completely explicit constructions tend to be computationally complex. We suspect that 
major progress on this point would come from developing other ways of understanding 
Kahler moduli stabilization. For example, mirror symmetry suggests the existence of type 
IIB non-Calabi-Yau deformations dual to turning on IIA NS-NS flux, which should thus 
be described by a similar superpotential depending on the IIB Kahler moduli. In this case, 
the Kahler sector would also contribute significantly to the counting of vacua, but since it 
would be governed by a flux type superpotential, it would also fit into the general class of 
supergravity ensembles for which the analysis below should be valid. We discuss this point 
a bit further in section 3.2.2. 

In any case, this discussion, together with the fact that IIB flux vacua are computa- 
tionally very accessible due to the underlying special geometry structure, justifies the claim 
that such ensembles are good models for the statistics of string vacua. 

One could do a similar computation for any class of models in which /C and W are 
computable. Generically, one would expect to obtain a similar distribution, in which a sum 
over K flux (or other) parameters leads to a rough independence for the first K quantities 
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in the distribution. In theories (such as the heterotic string) with fewer fluxes than the lib 
string, this would lead to constraints between the Z's and the (W, F) variables, while in 
theories with more parameters (such as F theory on fourfolds) one has more parameters 
and thus might expect a more generic ensemble of matrices Zjj. 

While it remains to be seen which distribution best represents the full set of string/M 
theory vacua, it seems plausible that this would be the one with the most free parameters, 
i.e. the fourfolds. One furthermore expects quantum corrections to these classical flux 
superpotentials. Thus another interesting ensemble, which might usefully represent this 
more generic distribution, is simply to take all of the parameters to be independent, 

= (fWcf"'FAd"''^"'+^^ZAB5{L - + |F|2 - |Zp) (2.12) 

where the delta function is an ansatz for a cutoff analogous to the one which made the 
number of lib flux vacua finite. This will turn out to be related to distributions previously 
considered in random matrix theory, as we discuss in section 2.4. 

To complete the specification of such a model ensemble, one must define the variables 
Uabc- Now in the fourfolds and other examples we might model this on, the U variables 
would still be determined in terms of the others, which might be important. One might also 
argue that a better model ensemble would be a simple proposal which keeps more of the 
structure of Calabi-Yau moduli spaces and their degenerations. These are interesting ideas 
for future work, but as we will argue below Eq. ( p. 12 ) already shares interesting features 
with more realistic distributions. 

In any case, we will keep most of the discussion general, and not assume any a priori 
relation between the variables in Eq. (|2.9|). 



2.2 Solving the equations V = Q 

The equations dV = are quadratic in the parameters, and the best way to exhibit their 
structure is to rewrite them as a matrix NiyV, Z) depending on (W, W ^ Z^ Z) acting on the 
vector Fa, as 



BaV I I Zab -2W J \f 



-2\W\ e^'^ZAB 
e'^ZAB -2\W\ 




(2.14) 



where e^^ = W/\W\. 

Define the matrix N(W, Z) to be the second of these matrices (i.e. the matrix depend- 
ing on \ W\), and write 

N = M - 2|I^| • 1 

in terms of the matrix 

\e''ZAB J ^ ' 
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Since Zj^b is symmetric, the matrix M is hermitian, so it has an orthonormal basis of 
eigenvectors. Thus, we see that non-supersymmetric vacua (solutions of = with 
-F 7^ 0) correspond to eigenvectors of the matrix M with eigenvalue 2|Vl^|. 

Clearly the matrix M is not a generic hermitian matrix; it has additional symme- 
try properties. First, its eigenvectors come in pairs with opposite eigenvalues (-|-Aa, — Aa). 
These eigenvalues are independent of Q\ this follows because their squares are the eigenval- 
ues of the hermitian matrix ZZ. The corresponding eigenvectors are 

where V'a solves 

Zi>a = Aa^a (2.17) 

and we take Aq > 0. Because the V'a are eigenvectors of the hermitian matrix ZZ, we can 
take them to be orthonormal. Then, defining the unitary matrix U by UAb = ipb,A, we have 

Z = UXU^ (2.18) 

where A = diag(A(j). Any symmetric complex matrix can be decomposed in this way. 
Thus, the generic super symmetry-breaking solutions to ( p. 13 ) can be written as 



2W = Xae''- F = fe^'/^iPa (2.19) 
with / G M. Varying 9 in [0, 2tt], this fills out a one complex dimensional subspace. 
2.3 Masses of moduli 

We have just seen that non-supersymmetric solutions of y = correspond to eigenvectors 
of a matrix M with eigenvalue 2|VF|. Thus, in ensembles such as the lib flux ensemble in 
which the parameters W and Fa can be freely varied, such vacua are generic. 

Of course, the vacua of most interest are metastable, i.e. the bosonic mass matrix V" 
has no negative eigenvalues. This mass matrix is 

(fV = (M + \W\){M - 2\W\) + Vl' + V2 (2.20) 

where the successive terms are of order F^, F^ and F'^. 

The F^ term is the product of the matrix N above with the matrix 

H = M + \W\-1 

which is the same as the matrix whose determinant enters into the density of super- 

symmetric vacua |11|. The higher order terms are 
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and 

^2" = ( ? ^ 1 . '^2 = RabcdF^F'' + Sab\F\^ - FaFb. 




We now assume \F\ << Mp = 1. In this case, the eigenvectors of M defined in 



( 2.16 ), are approximate eigenvectors of d?V. Let us take the eigenvalues Xa to be ordered, 

< Ai < A2 < ••• < Am- (2.21) 
The values of V" in the directions are 

= i±Xa + \W\){±\a - 2\W\) ± 2Re (e'^i^aSM + 2i^aS2i^a- (2.22) 



By (|2.13| ), one of these eigenvectors, call it ^p, is proportional to ^ e ^^F^ gidpB 
with A = 2\W\. There will also be a complementary eigenvector with eigenvalue — 2|1^|. 

Let us first consider variations of the moduli proportional to the other eigenvectors. 
Suppose the corresponding eigenvalue \a is greater than 2|1^|. In this case, the first term 
in ( ^.22 ) will be positive, and if the gap Aq — 2\W\ is sufficiently large (as will be the 



case for generic Z when F is small), this term dominates the higher order terms. On the 
other hand, positive eigenvalues less than 2\W\ will produce a negative m^, unless they 
are very close to 2\W\ and the higher order terms are fine tuned to compensate for the 
negative lowest order term. Thus, the bulk of tachyon-free nonsupersymmetric vacua with 
sufficiently small F will come from solutions for which 2|1^| equals the lowest eigenvalue 
Ai. 

We now consider (rra^)^ = (m^)^. For m^, the first term in Eq. ( p.22|) vanishes, so V" 
in this direction is given by the matrix element of the higher order terms: 

(m+)2 = ^ {Reie'^'UABcF^F^F'') + RaScdF^F^ F^ F^) . (2.23) 
Similarly, the matrix element for the complementary eigenvector with eigenvalue — 2|VF| is 
(nipf = 4\Wf + {-Re{e^''UABcF^F^F^) + RabcdF^F^ F^ F'') . (2.24) 

There are now two cases to distinguish. The simplest case is \F\ « \W\, but this can 
only be compatible with V = \F\'^ — 3\W\'^ + |Dp ~ if \D\ ~ \W\, i.e. the supersymmetry 
breaking is dominated by the D terms. If we have D = or even \D\ ~ \F\, then \F\ ~ \ W\ 
and we must take this into account. 

In the first case, and for generic Uabc, the first term in Eq. (2.23) will dominate, so 



(m^)^ will be positive for half the integration domain of 9. Furthermore, if \ W\ is not too 
small (greater than 0(\/F)), the first term in Eq. (2.24) will dominate and {rrip)'^ will also 



be positive. Hence the requirement of metastability puts only a mild constraint on the 
integration domain in this case. 



On the other hand, if \F\ ~ \W\, the first term in Eq. (|2.24 ) is of order F and no 



longer dominates. Instead, for generic U and 0, either {nip) < or (mp) < 0, as the 
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0{F) term appears with different signs in these. Thus, to get a metastable vacuum with 
zero (or positive) cosmological constant, we have to fine-tune the 0{F) term such that it 
becomes smaller than the 0{F'^) term: 

Ke{e'^/'^UABci^^^i^^) < 0{F). (2.25) 

This can be achieved by tuning 9 in an interval of size ~ which, if the measure for 
is otherwise uniform, will give an additional suppression of the expected number of vacua 
by a factor of \F\. 

However, because the diagonal mass matrix elements are now tuned one order in F 
smaller, we have to be a bit more careful and also consider the off diagonal matrix element 
between and ? This is 

{m^.cPV^t) = 2Im(e*^Vii5iVii) = -^\ui{e^'^UABcF^F^F^)- (2.26) 

Denoting si = e^^-i/^iSi-i/ii, ,82 = V'i/S'2^1, the determinant of the 2x2 mass matrix in the 
1-directions is thus 

det = 4 {\W\^{Resi + sa) - |si|2 + s^) . (2.27) 

Note that \i W ^ F, the leading term in the expansion in powers of F is generically 
— ~ -F^, and this is always negative independent of 6. Therefore, to avoid a tachyon, 
we really need < 0{F)'^, so (|2.25| ) gets replaced by 

\UABci^t^i^?\ < 0{F). (2.28) 

Assuming approximate uniform distribution of this complex component of f/, this metasta- 
bility constraint will therefore give an additional O(-F^) suppression rather than the 0{F) 
we deduced neglecting the off-diagonal element. The masses on the other hand will 
still be of order F. 

In any case, the overall conclusion of this analysis is that metastability is a rela- 
tively mild constraint, in the sense that it does not drastically reduce the number of non- 
supersymmetric vacua. In general, we find that a rough fraction 1/n of all nonsupersym- 
metric critical points are metastable vacua (since we can only use one of the n eigenvalues). 
This is much larger than the naive estimate of 2~-^"', obtained by taking the mf indepen- 
dent and symmetrically distributed around zero. This is consistent with the intuition that, 
for \F\ « \W\ « 1, the situation is similar to that for global supersymmetry, in which 
stability is automatic. But the arguments here apply far more generally than in this limit. 

Nevertheless, metastability can significantly influence the distribution: the distribu- 
tion of the lowest eigenvalue can be different from that of an arbitrary eigenvalue, and 
constraints such as Eq. ( 2.28| ) are important. 



■^We thank Michael Dine, Deva O'Neil and Zheng Sun for pointing this out to us. 
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2.4 Degeneracies and the dimension of the solution space 



We have just seen that for generic parameters W and Z, the space of possible supersym- 
metry breaking parameters F is one complex dimensional, the eigenvector of M with the 
lowest positive eigenvalue, multiplied by a general phase factor. 

However, at special loci in Z space, the matrix M might have /c-fold degenerate eigen- 
values. In this case, the vector F will vary in a A; real dimensional space, and varying 6 will 
fill out a k + 1 real dimensional subspace. This is the structure which might lead to power 
law growth of the number of vacua with the supersymmetry breaking scale, as suggested 
in |35, p^ ]. 

An explicit example in which this is realized is the "anti-supersymmetric branch" of 
lib flux vacua, with Z = W = or equivalently imaginary anti-self-dual flux. Now these 
vacua are physically not interesting, because they have cosmological constant at least of 
the order of the string scale , and moreover it can be shown that they always have a 
tachyonic mode. But the number of these vacua does grow as a high power of F, so we 
should not immediately dismiss the idea that some physically sensible subset of the vacua 
behaves in the same way. 

However, if we ask what we need to get degenerate eigenvalues in other circumstances, 
we find that they are non-generic, in the sense that one must tune more than k — 1 pa- 
rameters to get a fc-fold eigenvalue degeneracy. Because of this, it will turn out that these 
"higher branches" of non-super symmetric vacua are of lower dimension than the primary 
branch, and thus will not contribute to the overall volume estimate, and thus to the leading 
large L asymptotics for the number of vacua. 

The arguments are analogous to the more familiar discussion of degenerate eigenvalues 
in families of hermitian matrices, so let us review this first. For hermitian matrices, one 
must tune k'^ — 1 real parameters to get a fc-fold degeneracy. The simplest way to see this 
is to consider the change of variables to eigenvalues and eigenvectors M = U^XU. This is 
generically unambiguous up to permutation of eigenvalues, and up to a {/(I)" left action 
on U. Thus the real parameters of M go over to n eigenvalues and — n coordinates on 
U{n)/U{l)"' . However, a diagonal matrix with a /c-fold degenerate eigenvalue is preserved 
by conjugation by an element of U(k), and thus such matrices form a stratum in the larger 
space of all matrices of real codimension k — 1 (required to tune the eigenvalues) plus k"^ — k 
(the dimension of the stabilizer group). 

To make the analogous argument for the case at hand, we consider the change of 
variables Eq. ( 2.1^ ). Now the n(n + 1) real parameters in Z are rewritten as n eigenvalues 
and V? parameters of the U{n) group element V . On the other hand, a A;-fold degenerate 
diagonal matrix will be preserved by a unitary V satisfying VV''" = 1, in other words an 
element of SO{k). Thus the degenerate matrices form a stratum of total real codimension 
{k + 2)[k — l)/2, the dimension k{k — l)/2 of SO{k) plus k — 1 for the eigenvalues. 

A related way to see this is to consider the change of variables from the Lebesgue 
measure on matrix elements, to eigenvalues and unitary group elements. For the hermitian 
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matrix, this is 

c/^'m = [dU] J] d\i - \jf 

i i<j 

where [dU] is the invariant measure on U{n). This exhibits the extra real codimension 
k{k — 1) as an exphcit scahng dimension in the Jacobian, in analogy to the scaling r^~^dr 
of Lebesgue measure in polar coordinates. 

The analogous matrix ensemble for the present case is determined by the symmetry 
properties of Eq. ( |2.15D . Introducing a triplet of Pauli matrices dj to act on the 2x2 block 



structure of Eq. ( ^.15 ), these are 

M = M^ = -asMas = -a2M*a2. (2.29) 

These symmetry properties are formally equivalent to those of the "CI ensemble" of Altland 
and Zirnbauer |j^] . This was introduced as an ensemble of Hamiltonians of mesoscopic sys- 
tems, and in that context these symmetry properties are time reversal invariance and spin 
rotation invariance. While the physics and the interpretation of the symmetry properties is 
different here, and the actual ensembles arising from string theory (as discussed in section 
2.1) are a subset of the CI ensemble, certain properties of the CI ensemble and specifically 
the structure of eigenvalue degeneracies should be shared by the stringy ensembles. 

Thus, we consider the CI ensemble measure, in which the different matrix elements 
(consistent with Eq. ( p. 29 )) are independently distributed; this is just the ensemble Eq. (|2.12| ) 



We then perform the change of variables 

M = 




The CI measure transforms to 

dfi[\, U] = [dU] n d{Xl) H \Xl - Xl\ (2.30) 



a<b 



and we again exhibit the codimension k{k — l)/2 + (k — 1) of the symmetry locus as a 
scaling dimension, as well as the fact that the true invariants of M are the squares of the 
eigenvalues. 

Now, given a set of k eigenvectors ipa, 1 < a < A;, all with eigenvalue A, the analog of 



Eq. (|2T9D is 

2W = Ae^^ F = Te^^/^a (2.31) 

where f"" is a real k component vector. Thus, while the associated locus of nonsupersym- 
metric vacua is + 1 real dimensional, and the integral over fluxes d'^^F would lead to 
an additional power law growth iFj'^"^ (compared to the case A; = 1), the variables Zjj 
must satisfy (A; + 2)(A; — l)/2 additional real constraints, leading to a total real codimension 
{k + 2){k — l)/2 — (fc — 1) = k{k — l)/2 for the fc-fold degenerate branch. Thus already for 
k = 2 these branches are of lower dimension in the generic case, and will have volume zero. 

A similar discussion applies in the special case of a A;-fold degenerate zero eigenvalue, 
but now both ^+ and ^~ can be superposed in Eq. ( |2.19| ), so the resulting locus of 
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nonsupersymmetric vacua is 2k real dimensional. On the other hand, a /c-fold degeneracy 
at zero generically appears in real codimension 2{k — 1) + k{k — 1) = (A; + 2)(k — 1), so 
again the resulting branch is of lower dimension. 

There is a loophole in these arguments, which is the assumption of genericity. If 
the matrix Zjj depends on the fluxes and fields in a non-generic way, such that a A;-fold 
degeneracy appears with real codimension k — 1, then the corresponding branch of vacua 
would have the same dimension as the primary branch. This is what happens in the "anti- 
supersymmetric" branch: because of the relation Zjj = J^jjkZ^ , one need tune only the 
n complex parameters Zqj to obtain Z = and a 2n-fold degeneracy at zero, so this branch 
is again of the full dimension. Since the matrices Zjj obtained in lib flux compactification 
are definitely not generic, we need to understand this point. 

A familiar way to get a k-fold eigenvalue degeneracy in codimension k — 1 for a family 
of hermitian matrices is to take the sum of k commuting matrices, 

M = tiMi with [Mi, Mj] = 0. (2.32) 

In this case, since can simultaneously diagonalize the Mi, the problem reduces to con- 
sidering families of eigenvalues, for which there is no difficulty in tuning to degeneracies. 
Conversely, if the family is a linear sum Eq. ( p. 32 ), this is the only way to get eigenvalue 



degeneracies; if some [Mj, Mj] ^ then eigenvalue repulsion always makes the codimension 
of a fc-fold degeneracy (reachable by tuning parameters which include t j and tj ) higher than 
A; - 1.3 

The situation for matrices in the symmetry class CI is similar, and governed by the 
same criterion Eq. ( p.32| ) applied to the matrix Eq. ( 2.15| ). This implies that a fe-fold 



eigenvalue degeneracy in codimension k — 1. is possible only if Z is a linear sum of k 
matrices Zi which commute after applying a similarity transformation; in other words if 
there exists a unitary U such that 

[UZ^U^,UZjU^]=0. (2.33) 

Thus generic sums of n matrices will exhibit the generic behavior of an n-parameter family 
of matrices, in other words eigenvalue repulsion and no n + 1-fold degenerate eigenvalues. 
The most obvious way to get the structure Eq. ( p. 32 ) and independent F breaking 



parameters, is the case that the effective supergravity Lagrangian can be written as the 
sum of two (or more) completely independent Lagrangians, each depending on a distinct 
subset of the fields. This fits with the intuition that in a model with several decoupled 
hidden sectors, supersymmetry breaking in the various sectors would be independent, but 
only if the different sectors satisfy a very strong decoupling requirement. One way to get 
this is to consider limits in moduli space, and following up this idea leads to the question, 
what would be the consequences of a nearly degenerate eigenvalue. From Eq. ( |2.20D , this 



^This can be seen by considering the free classical mechanics with Lagrangian tr (dM/dt)^, and the 
solution M = (1 — t)Afi + tA'h. Changing variables to eigenvalues and unitaries, if [A/i, M2] 7^ there will 
be a non-zero angular momentum L — [M, M], and a potential Lfj — repelling the trajectory from 
the degeneracy. 
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would seem to lead to a nearly massless boson, but no obvious enhancement of the number 
of vacua. This case deserves closer study, but we leave this for future work. 

It turns out that less simple but naturally occuring examples can also satisfy Eq. ( 2.32| ) 
and lead to degenerate eigenvalues; we will see this below for the case T^/Z2 x Z2. However, 
this example has an unusual degree of symmetry, and we have seen no sign of the structure 
Eq. (2.32) in other Calabi-Yau compactifications. 



2.5 Some physical comments 

While we suspect that some of these observations have been made before, our discussion 
does not look much like the existing discussions of supersymmetry breaking in the literature, 
which made it hard for us to compare with previous work. There are two reasons for 
this. First, in model building terms, the considerations here would describe the origin 
of supersymmetry breaking in hidden sectors, and not its mediation or effects on the 
observable sector. Second, we have tried hard to use the underlying geometry to simplify 
the problem to its essentials, perhaps at the cost of some physical intuition. Let us make 
a few comments to remedy this. 

The matrix M is closely related to the supergravity fermion mass matrix, which can 
be read off from the Lagrangian ([p^, (23.3)): 



where is the gravitino, XA are the fermionic partners of the moduli , and i/^"" and XA 
are related to these by complex conjugation. Making the change of variables g«6»/2^a 
and XA e~*^''^XA; and decomposing x into the eigenvectors Eq. ( 2.17 ), one finds that the 



eigenvector ^/JF with eigenvalue 2\W\ is the goldstino, while the other eigenvectors "i/^a lead 
to fermions with Major ana masses niA = Xa- The condition that M has an eigenvalue 2\W\ 
is also the same as the condition for a scalar in a supersymmetric vacuum to go massless, 



as discussed in section 3.2 of |11]. This is not a coincidence; it is because different physical 
branches of the space of vacua attach at second order phase transitions, i.e. points at 
which a scalar field becomes massless. As we vary A^, when its kernel jumps, a modulus 
will go massless. Note however that there are two types of variation of N, those which 
come from varying the moduli z^, but also those which come from varying the effective 
Lagrangian data {W,F,Z). While all supersymmetry breaking branches are connected to 
the supersymmetric branches by varying A^, these will be physical second order transitions 
only if they are connected by varying moduli. 

One way to see the origin of the stability condition > Ai = 2\W\ is that the 
bosonic mass matrix Eq. ( |2.6D has a universal term dy^^di^'p — 2|VFp), which for |F| << 
\W\ is a universal tachyonic contribution. This is not much emphasized in the literature, 
which usually considers the case D = or equivalently |Fp = 3|VFp, in which case this 
contribution is positive. However it is quite important in general, especially in the case 
|F| « \D\ ~ \W\. 
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3. Distributions of non-supersymmetric vacua 

We proceed to work out the distribution of supersymmetry breaking parameters, 

d^[F] = j dfi[W,F,Z,U] 6'^"'{dV)\detd^V\. (3.1) 

We start by changing variables from Zab to its eigenvalues and eigenvectors Xa and tpA- 
This gives us a joint distribution 

dfi[W,F,XA,iPA,U]. 

These variables are redundant under permutations in Sn acting on the indices A. To fix 
this, we take the Xa in increasing order, as in Eq. ( 2.21 ). 

We then solve the constraint dV = as in ( p. 19 ), in a basis with Z diagonal. Put 
(l)A = arg(e~*^/^i<A)- We focus on solutions with 2\W\ = Ai, because this is where the bulk 
of the metastable vacua are located, as argued in section 2.3. Factorizing 6{N ■ F) (with 
N -F asm Eq. ( pT3D ) turns 

^2mp 5(2m)(^ . ^ ^2m p J^s[Re{XAFA - 2\W\e-^^ Fa)] 5[Iui{XaFa - 2\W\e-^^FA)] 

A 

= d^^'F 1[5[\Fa\ cos <PA{XA-2\Wm\FA\ sin ^a{Xa + 2\W\)] 
A 

_ df 5{Xl-m') ,,,, 
" I/I Idet'iVl ^'-'^ 

where det'A^ is the determinant with the modes ^'^ excluded, and F is given by ( 2.ig| ) 
with A = l. 

For small F, the determinant of di^V can be approximated by 

detdV ^ X\{m\f{m-^f 

A 

« (m+)2(mr)^ \[i\W?- >^\mW? - Ai) 

■det'iV(m+)2(m^)^ (3.3) 



A>1 

detH 
'3\WV 



The niA are the masses introduced in Eq. ( p. 22 ), and we assumed the Xa for ^ > 1 to be 
sufficiently above Ai = 2|VF| such that the higher order terms can be neglected in these niA- 
Thus we can extract the det'A^ from the numerator, to cancel with the one in Eq. (p.2|). 
The l/3|VFp here cancels the factor in detH coming from the eigenvectors so this 
expression is non-singular. 

The result for the distribution of nonsupersymmetric vacua is dfj,[W, F, X,ip,U] p, where 

p = 5'--\F) 6{Xl - m?) • ^^^^l^p^ (3.4) 
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and 5^"^ ^(F) is interpreted as above, and {mf)'^ is as in Eq. (|]2|) and Eq. d^H) 



Most of the structure of the result is in the factor | det H\. This is the same as the Ja- 
cobian appearing in the density of super symmetric vacua, so the result can be summarized 
as saying that the density of non-supersymmetric vacua is very similar to that for super- 
symmetric vacua, with a "correction factor" proportional to the product of the squared 
masses of the two moduli in the supersymmetry breaking direction. 

If we were on a locus with k degenerate eigenvalues, we would have to generalize the 
above steps: Eq. (^) would contain a product of k delta functions 5(Aj — 2|VF|)/F, we 
would have a larger kernel of N, and would need to replace Eq. (^^) with the determinant 
oi a k X k submatrix of V^' + V2 ■ While possible in principle, normally this will happen 
only in codimension larger than k, as discussed earlier, and thus will not contribute to the 
final integral. 

If we do not impose conditions on the value of the cosmological constant, we have 
generically (i.e. for \W\ not too small) (m-^)^ ~ 4|VFp and the distribution becomes 

Here 0+(^) restricts 9 to values for which [mf)'^ > (which is approximately half its 
integration domain). Thus we see that the small \F\ behavior of the distribution is dF, 
so the number of generic metastable supersymmetry breaking vacua with supersymmetry 
breaking scale Eq. ( |1.1D less than M^: 

Af{Msusy < M,) ~ M2 (3.6) 

If on the other hand we restrict to metastable vacua with 1/ ~ (and D = 0),^ we 
have to multiply this distribution by 6{\F\'^ - 3\W\'^) = 6{\F\ - \/3\W\)/2\F\, leading to 

p = 6{\F\ -V3\W\) 5^--\F) 6{Xl - A\W\') 0^(0, F) • ^"^^J^J^^^' (3.7) 

where 0+(^, F) restricts 9 and F to values such that (mf )^ > 0, which as argued in section 
2.3 gives at least an additional 0{F) suppression. We also used | det H\ « 3|M^p nyi>i ^'a ~ 



f |detZ|2. 

Let us now grant that the joint distribution we started with looks like 

dp[W, F, Z, U] ~ d^W d^'^^^F dfi...; 

in other words W and each component F^ are roughly uniformly distributed in the complex 
plane, and independent of each other and the Z's. All this is true in the lib case. We can 
then use the delta functions to solve for W and \F\, to obtain 

p~dMA.C/]dF9,(<».F)l^. W';^"-i')' . (3.8) 



*By ~ we mean \V\ « 0{F), which can stiU be much larger than the observed cosmological 
constant. In particular further quantum corrections to A remain approximately within this window. 



- 15 - 



As discussed in section 2^, metastability forces (m^)^ ~ {m^ )^ ~ F'^ and Q-^-{6,F) ~ F^, 
leading to 

p ~ dplXlF'^dF. 

Furthermore, in a generic distribution of superpotentials such as Eq. ( p.30 ), one expects 

dfi[X] ~ XdX ~ FdF, 

leading to 

p ~ F^dF. 

Thus, the number of metastable susy breaking vacua with near-zero cosmological con- 
stant and susy breaking scale Mgusy < goes like 

M{M,usy < M* , A ~ 0) ~ (3.9) 

Among such vacua, small supersymmetry breaking scales are disfavored. Note also that 
this result suggests that metastable de Sitter vacua obtained by pure F-term susy breaking 
are relatively rare. 

3.1 Type IIB flux vacua 

For type IIB flux vacua (ignoring Kahler moduli), we have [11 1: 

Zoo = 0) Zoi = Zj, Zjj = J^ijkZ^ (3.10) 

and 

Uooi = 0, Uoij = J'ijkF^, Uijk = DiJ^jklZ^ + J'ijkF^. (3.11) 
The nonzero curvature components are 

-^0000 = ~2, RlJKL = ^IKM^^ JL — ^IJ^KL — SuSjK- (3-12) 

This gives 

{mff = (\ruKF'F^ ± 2e-^^'FoFj\^ - 2\F\^ ± Re {e^'' DjJ^jklZ' F' F^ F^)) 

+5±,_4|iy|2, (3.13) 

so indeed when F ^ W and DT ~ 0(1), to keep both squared masses positive, one 
needs an order F fine-tuning of the term proportional to DJ^. At large complex structure, 
DJ^ = so no such fine tuning is needed.^ On the other hand, if the number of moduli 
is large, this constitutes only a tiny part of the moduli space and hence of the number of 
vacua. Near a conifold degeneration, both DJ^ and J- blow up, with DJ^ ~ J^'^, so the 
same kind of fine-tuning is needed as in the generic case. 



'However, the remaining terms may still fail to be positive. 
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3.1.1 One complex structure modulus 

We turn to some more explicit results, beginning with the simplest case of one complex 
structure modulus. We will see that this illustrates some but not all features of the generic 
discussion above. 

Letting Z = Zqi, the matrix Zjj = DjDjW is 

Zij = 

and the eigenvalues A of M satisfy 

(A^ - |Z|2)(a2 - |Z|2(1 + |^|2)) - = 




More explicitly, Ai^2 = ^{\^\ =t + 4)|Z|. Near a conifold limit JF diverges, and 

Ai ~ IZI/I^I, A2 ~ \J^\\Z\. The measure for Z in this limit is d\Z\'^ = jJ^pdAf. 

As in the general discussion, we consider |T^|,|F| << 1. Here this requires taking 
|:F| >> 1. The constraint in Eq. (|2.1l|) forces |Zp ~ L, so we write Z = -/Le*''^. Then, at 



a metastable pure F-breaking critical point of V, |-F|V3 ~ = Af/4 ~ L/4\T\'^, which 
is indeed small when !F ^ oo. Granting ~ F"^ and substituting into Eq. ( |3.8[ ), we find 



d^l[F] =jd'w d'z dF de 6i\z\^-m\F\ - V3\w\)6{xl-i\w\')e+ie,F) ^ ■ ^^^^^^ 

~ d(j) de Q+{e, F) X dF 5{F - 

Granting the behavior G+ ~ \F\'^, this looks like a scaling F^dF at a given point in moduli 
space. What happened to the \F\^dF we found earlier? One factor was cancelled by 
the in the denominator - the general expectation that | detZp would go as \W\'^ is 
violated in this case, because the see-saw mechanism pairs a large eigenvalue with the small 
eigenvalue. The remaining missing factor of F comes from the fact that in this particular 
ensemble, there is no additional suppression of small A from the Z-measure, as Z does need 
to be tuned small to make A small. 

Finally, we need to incorporate the dependence of J- on moduli. Suppose we are near 
a conifold limit, parameterized by a single complex structure modulus t — > 0; then 

~ ~ ~ LlJ'l'^ ~ LIA"^ log^ (3.14) 
and the measure on moduli space d^zy/detg becomes 



d^t log lip ~ 



Llog^ |F|2 
leading to 

The final distribution (with the metastability constraint leading to 0+ ~ F"^) goes as 
F^dF. The difference with our previous general claim F^dF arises from the see-saw factor 
1/1 in I det (required by the scaling det Z ~ F^). 
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Thus, despite the fact that flux vacua are dual to gauge theory in the conifold limit, 
and do lead to hierarchically small scales, the final vacuum distribution does not show 
a corresponding enhancement at small scales. This came from a combination of effects, 
which we summarize again here and in the conclusions. 

The hierarchically small scale is \t\, and it is true (Eq. ( |3.14| )) that this leads toW^t, 
which after imposing A = implies F ~ t. On the other hand, the enhancement of 
the number of vacua near the conifold point found in [^] is not present for the small W 
component, in present terms because the large matrix element cancels out of det Z. In 
addition, the metastability condition Eq. ( 2.28| ) gives an extra F^. 



3.1.2 A simple but nongeneric multiparameter example: T^/Z^ 

Another example which can be treated exactly is the orbifold {T'^)'^/Z,2 (without discrete 
torsion). This has three complex structure moduli Tj, and the only nonzero J-ijk is .^^123 = 
1. As usual we ignore the Kahler moduli. The Z-matrix is 



Z 



( Zi Z2 Z3 \ 

Z\ Z3 Z2 

Z2 Z^ Zi 

V Z3 Z2 Zi y 



Its eigenvalues and eigenvectors can be obtained by elementary means. The eigenvalues 
are 



Ai — IZ1I+IZ2I — IZ3I, A2 — |Zi| — IZ2I + IZ3I, A3 — — IZ1I + IZ2I+IZ3I, A4 — IZ1I+IZ2I + IZ3I, 
and the eigenvectors tpA have simple expressions which depend only on the phases (pi of Z^: 



(l -1-1 1 \ 



I- 11-1 

II- 1-1 
\i i i i 



\ 











3i{9il-02-<^3)/2 











\ 



Q g«{~<^l-<^2+03)/2 



Note that since the Z-measure is just d^Z, there is no suppression of coincident or zero 



eigenvalues. This can be understood along the lines of the discussion at the end of 2.4, by 
noting that the matrices Zj = dZ/d\Zi\ satisfy Eq. (|2.33| ) with = {iPa)^ independent 
of |Zj|, and thus can be simultaneously diagonalized on a three real dimensional slice of 
parameter space. 

This feature will be destroyed by generic perturbations of the prepotential. A simple 
example for which this is easily verified is obtained by adding "conifold-like" terms J-ijk = 
SijSiK/tR- Thus, we believe the present example is not typical, and that more generic 
flux ensembles will share the CI ensemble behavior. However, let us finish the discussion 
for completeness. 



Putting F = ipi, Eq. (|3.13D gives 



2\F\ 



im^y = 4i\W\' 
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This simplification and in particular independence of 9 and vanishing of the 0(F) term 
occurs because DJ^ = in this case. This is another nongeneric feature of this model. 
One consequence is that there is less tuning control of the masses, and in particular in 
the case of pure F-breaking and approximately vanishing cosmological constant, (rn.^)^ = 
-8|F|V3 < 0, so there are no such metastable vacua. When a constant D-term is added, 
this becomes {mi)'^ = 4{D'^ — 2\F\'^)/3, which is positive for sufficiently small F. 

Plugging this in the general formulae, we find for the vacuum density at zero cosmo- 
logical constant 



dfi ~ dAd^W d'^Z dF 5{\F\^ + - 3\W\^) 6{Xi - i\W\ 



xe{\w\ - \F\) . {x,x-A,f . ^-^^^^ 



^dAdF\-W=MlL^^\ 

The exponent of L^: can be checked by counting dimensions, where F, W, Z get assigned 
dimension 1, and thus L, A get dimension 2. The total dimension must be 263 = 16, which 
is indeed the case with the above power of L^:. 

When F « D ^ W, this goes as dFF, which is as in the generic case but with an 
extra factor of F, due to the fact that (m^)^ is of order F"^ here. When F D (but not 
too close to D/\/2), this goes as dFF"^. This has two powers less than the generic case, one 
because there is no additional suppression here of small A, and one because no fine-tuning 
is needed to kill the ±0{F) term in the masses. Finally, when \F\ —>■ D/y/2, the density 
drops to zero together with the mass . 

Finally, there would also be terms from the branches with degenerate eigenvalues, 
of the same form multiplied by powers of Fl\fTj. If we trust the results all the way to 
|F| ~ \fTjMl^^^^^^ which might be reasonable for Mstring << Mp (as at very weak string 
coupling), these will be comparable in number to the vacua we just described; however as 
we discussed their existence appears to be a special feature of this example. 



3.1.3 The generic multiparameter case 

We now discuss what we would expect to see for a more generic multiparameter case, in 
the continuous flux approximation. To get low scale supersymmetry breaking with zero 
c.c, we need small W amd F. While in this approximation, there is no a priori constraint 
on the size of these parameters, solving the equations Eq. (|2.13| ) requires the matrix Z to 
have a small eigenvalue. 

It is hard to tune a generic matrix to get a small eigenvalue. This is easily made 
precise in the standard ensembles of random matrices. In the CI ensemble, the measure 
for matrices with a small eigenvalue A goes as (i(A^). Following this through leads to the 
generic F^dF distribution we discussed above. 

In the one parameter case, Z is a 2 x 2 matrix, and the dilaton dependence of the 
lib flux superpotential leads to a see-saw structure. This is special to 2 x 2 matrices. 
More generally, what we can accomplish by tuning to conifold points is to get large matrix 
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elements in Z, for example the TZ ~ Zjt. In higher dimensions, large matrix elements do 
not generally lead to small eigenvalues. 

A simple model which may illustrate the general situation is to take Z to have generic 
order 1 coefficients except for a single large matrix element Z\\ = T » 1. The eigenvalues 
of such a matrix can be found by treating the 0(1) coefficients as a perturbation around the 
spectrum of the matrix for which only Z\\ = T \s non-zero. A generic 0(1) perturbation 
will shift these eigenvalues by 0(1), and the resulting eigenvalues will be .7-" + 0(1) and the 
rest 0(1). 

The see-saw matrix escapes this general result because Zqo = forces det Z ~ 1 and 
thus a small eigenvalue must appear. To get this effect in a higher dimensional matrix, 
the determinant of the 1, 1 minor of Z must vanish. This would appear to be a rather 
complicated and non-generic tuning. 

Thus, we see no clear way out of the generic F^dF prediction of our earlier discussion 
in this case. Since even in the one parameter case, the expected enhancement of low scale 
vacua was cancelled by measure factors, the suppression of low scale vacua of this type 
would appear quite general. 

3.2 Quantized fluxes 

As in IQ, |ll[, we have ignored flux quantization in the results so far, instead taking the 
fluxes to be continuous variables. At first sight this may seem a drastic simplification as the 
allowed values of the EFT parameters {W, F, Z) are heavily influenced by flux quantization. 
For example, the hierarchically small value of certain contributions to the superpotential 
near the conifold point is only apparent in the quantized flux problem. 

Now, as in |Q, |ll|], one can argue on general grounds that the approximation of taking 
flux continuous reproduces the leading large L asymptotics for the number of vacua, because 
the volume of the region in flux space supporting vacua is the asymptotic for the number 
of quantized flux vacua (lattice points contained in the region). We make this argument 



below. Furthermore, explicit numerical study |27| , 11, 10, |1^ has confirmed the validity 
of the approximation in counting supersymmetric vacua when L > K/r'^, to estimate the 
number of vacua in a region of moduli space of radius r, in finding the distribution of W, 
and other observables. 

This does not mean that we should immediately accept the analogous claims for non- 
supersymmetric vacua, as there are clearly important differences between the two prob- 
lems. For example, the validity of statistical approximations to obtain the W distribution 
is surely helped by the fact that is a sum of complex quantities with arbitrary phases, 
which typically involves many cancellations. On the other hand, since the quantity M^^^y 
is a sum of squares, errors in the approximation might lead to systematic overestimates. 

More work will be needed to find the minimal fiux L for which the results here are 
accurate, but let us proceed to make some general comments. 

3.2.1 Geometry in flux space 

We now describe the region in type lib fiux space which supports nonsupersymmetric 
vacua, in the same sense developed for supersymmetric vacua in 0. 
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We work around a point z in moduli space and make the appropriate decomposition 
iyV, -F, Z) of the fluxes at that point. We then think of the (W, F) directions as fibered 
over the Z plane. At a given Z and 0, the branches of nonsupersymmetric vacua sit at the 
m values of \W\ given by the eigenvalues of Z, and in one-dimensional subspaces of the F 
space given by the eigenvectors. These m solutions of the eigenvalue equation or sheets are 
the basic "branches of solutions" in this problem. Of these, only the one with the lowest 
value of \W\ is tachyon-free. Moving around in the Z plane varies them, and one can have 
monodromies about points with degenerate eigenvalues which exchange sheets. 

Varying then fills out a correlated circle in W and F, to give two real dimensional 
sheets in the 2m + 2 dimensional (VF, i^) space. Finally, varying the 2m moduli will turn 
these sheets into cones or balls of full dimension. 

Another way to think about this is to consider the "universal" solution space as (W, F) 
fibered over the entire m(m-|-l) real dimensional space of Z, which contains a real codimen- 
sion 2m domain of allowed nonsupersymmetric vacua (say with and jl^p bounded) 
which looks like a two dimensional sheet in {W^F). Each point in moduli space then pro- 
duces a rank 4m lattice (in the lib theory; on fourfolds it could be larger) in this space, 
and then varying the 2m real moduli allows these to hit the universal solution space at 
isolated points, the non-super symmetric vacua. 

In any case, the region in flux space which supports nonsupersymmetric vacua is of 
full dimension in the flux space and has a smooth boundary, and thus standard arguments 
imply that the leading large L asymptotic for the number of lattice points contained in 
this region will be its volume. 

3.2.2 Subleading components in L 

Perhaps a subleading term at large L has a different F distribution, for example producing 
many more small F vacua, in a way which makes it dominate in the physical regime. 

While the true distribution of F is quantized, one would expect this to cut out small 
F vacua. In particular, one might naively expect the distribution at a value F to be well 
approximated by the large L asymptotic only when L > On the other hand, small 

parameters enter into this relation, as we saw in the explicit one parameter example, so 
the situation in cases of interest is probably better than this. In any case, this effect will 
shift the F distribution towards higher scales. 



On the other hand, it has been suggested by Dine et al [15| that the W distribution 
could obtain a component highly peaked at zero, which after tuning the c.c. would lead to 
a peak in F. Their idea is that, because W = restores R symmetry, there might be an 
enhanced number of flux vacua with W = 0. Then, since our exact considerations are of 
course just approximations to the full physical problem, one might expect such a peak to 
be smoothed out to an enhancement of small \W\ vacua in the full theory. 



Recently DeWolfe et al |14] have studied the problem of finding super symmetric flux 
vacua with 1^ = in some detail and indeed find enhanced numbers of = vacua at 
subleading order in L in simple examples. While it remains to be seen how important this 
effect is, the idea and evidence for it seem quite reasonable. However, it seems unlikely 
that this comes with an enhancement of small W vacua in the pure flux vacuum problem. 
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This is for both mathematical reasons (the points with an enhanced number of vacua tend 
to be surrounded by "voids" without vacua), and physical reasons - while the flux vacua 
do contain exponentially small effects, the superpotentials which stabilize all moduli also 
contain other 0(1) contributions which show up in W and cannot be eliminated. 

Rather, one needs to call on additional physical effects to smooth out the distribution. 
At first sight the suggestion of |jl5[ seems plausible: exponentially small corrections could 
lead to vacua with small W ~ exp —1/g'^, and a uniform distribution of couplings dg would 
translate into a logarithmic distribution dW/W. However, at present we know of no explicit 
ensemble of models in which this idea would be realized. In particular, there is no evidence 



that the KKLT construction |33| would do this. While it relies on exponentially small effects 
to stabilize Kahler moduli, these are balanced against a preexisting small W obtained from 
the flux sector, and the construction does not work without this. One might hope that a 
model more like the original "racetrack," in which several competing exponentials stabilize 
the Kahler moduli, could lead to small W, but if the potential is naturally a polynomial 
in the exponentials qi = ex.p—1/gf, one might expect the resulting distributions to be 
uniformly distributed in the variables qi. All this is not to say the suggestion is clearly 
false, but rather that to support it one needs to show that it is realized in some explicit 
ensemble of models which could plausibly come out of string compactification. 

In our opinion, at present the best motivated conjecture one could make for such an 
ensemble of Kahler stabilized models, is simply that it is similar to one of the known flux 
ensembles. As in [|l8|, one might try to argue this from the existence of gauge- flux dualities 
such as [^] which (in much simpler examples) explicitly relate the two classes of vacua. 
Not having a compelling argument of this type, we will simply make the comment that if 
the true ensemble of Kahler stabilized lib models turns out to be different, one will also 
want to understand why this type of duality argument fails. 

Suppose we were to grant that there is a large component of vacua with small W 
distributed as d'^W/\W\'^; how would this influence our results? We need to first ask if 
this component is already visible in our computations, and properly taken into account. 
Although there are similar looking factors in our intermediate steps, in fact they are not 
present in the final results, for the reasons we explained. Rather, we would interpret the 



suggestion of |15] as saying that considerations in a sector of the theory not considered here 
lead to an additional factor |VF|~^ in the original distribution Eq. ( 2.11 ). If we incorporate 



such a factor, given F ^--^ W we will find the generic distribution changes from F^dF to 
F^dF. This stih appears to favor the high scale, at least for the purely F breaking vacua. 

4. Including D-terms 

At present there is no explicit ensemble of EFT's with D terms which looks simple and 
universal enough to justify a detailed study. Thus in this section we simply add a generic 
additional contribution to the potential, 

V= |F|2-3|T^|2 + eVD, 

and discuss its consequences. 
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To get small supersymmetry breaking, we take e to be small. This would be the case 
for instance for a contribution arising from an anti-D3 brane placed at the bottom of the 



warped throat developing near a conifold degeneration |2G]. In this case, which will be our 
concrete example below, e ~ where v is the complex structure coordinate given by 

the period of the vanishing 3-cycle. ^ 

The condition for a critical point W = becomes 




V{eV^) = -ed 



with as in Eq. ( 2.13| ). This is solved by 



eN-'d. (4.1) 



Decomposing d = d^^^ + d^'^^ with G M and as defined in Eq. ( p. 16 ), we have 

iXA-2\W\y ^ (Aa + 2|M^|) 



i7|2 = ,2V- (4)! , KJ! (4 2) 

I a . - 9iwh2 ^ a . -u 9iT/r/h2- 



A 

4.1 Lifted susy vacua 

Generically (i.e. for not too close to 2|H^|), we can assume all terms in the sum to 
be at most of order 1. This assumption needs some discussion for the antibrane example 
because of the dependence of e on the modulus v, which implies d^ ~ 1/v. However, the 
matrix element Z^^ ~ T^^^Z"" ~ l/u in compensates for this. In other words, there will 
be a pair of eigenvectors ^'^ approximately associated to the u-direction, and although in 
this direction d^ ~ 1/f, we also have ~ 1/u, and the two cancel against each other in 
Eq. ( |4.2| ).'' Thus, \F\ is generically of order e. 

In this case, the D-term dominates over the |Fp-term in the potential: 

V ^ -?>\W\^ ^eVB, 

and the measure becomes 

b''^(dV)\ det V"\ = '^""(^^ Idet'Ir^l^''"^^^^ ' '^^^ ^^^"^ + {HN)-\Vl' + + eVl^))\ 
« 6^"^ {Fa + ee'^{N-^d) a) \ det H\. (4.3) 

Again some discussion is needed to justify dropping the V^' term, since it blows up as 

in the u-direction. However in this direction HN ~ A^ ~ 1/u^ as well, so this cancels out. 



^While it has been argued ^ that if the anti-D3 contribution has an EFT description, this must be 
a D term, we know of no description within the usual rules of TV = 1 supergravity which satisfies all known 
properties of the string theory construction. However this point will not be crucial for what follows. 

^We should actually go to an orthonormal frame first, but since the metric ds^ ~ log \v\^^dvdv in the v 
direction, this only induces logarithmic corrections. 
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The above expression is exactly the density for the supersymmetric branch. Similarly, 
the condition for the absence of negative modes of V" is the same as for the superymmetric 
case: > 2\W\, up to possible corrections from V^, but if — 2\W\ » e these 
corrections will be small. 

Thus, these nonsupersymmetric vacua correspond to supersymmetric vacua "lifted" 
by the D-term. The number density of such vacua at cosmological constant A will ap- 
proximately be equal to the number density of supersymmetric vacua at susy cosmological 
constant A — eVn- Since the latter is nonvanishing at zero we have that for small 
A < eVo, this density is essentially independent of A, and roughly equal to l/MtotiT^) 
when integrated over a region IZ in moduli space. 

The supersymmetry breaking scale Mgusy for these vacua will be of order y^, and thus 
the expected number of these "lifted susy" vacua with Mgusy < Af* is, using the results of 



1 11] for the susy vacuum distribution near the conifold. 



M{Msusy < M,) ~ (4.4) 

That is, taking into account tuning of the Higgs mass, for this family of near-conifold lifted 
susy vacua, a low susy breaking scale is favored. 

4.2 Perturbed F-term susy breaking vacua 

The situation changes when one of the eigenvalues Xa approaches 2|1^|. Since we turned on 
the D-term potential as a small perturbation, we still expect that 2\W\ will be the generic 
approximate lower bound on the Xa_ to get a metastable minimum. Let us therefore assume 
that all Xa are well above 2\W\ except possibly Ai. Define u = Xi — 2\W\. When u <C 1, 
we have 

F^l^^l, (4.5) 

u u 

Note that since we are interested in <C 1, we require u ^ e. All (m^)'^ for j4 > 1 as 
defined in Eq. ( p. 22 ) will still be given approximately by their susy values. On the other 
hand 

(m+f = {u + 3\W\)u + ci- + + c+e (4.6) 

e 

{m^f = {u + \W\){u + 4\W\) -CI- + + cg e, (4.7) 

where the q are generically of order 1. We consider three regimes: 

(1) < u < 1: 

In this case, the first term in Eq. ( [4.6| ) is much bigger than e^^^, and the remaining 
terms are all much less than e^/^. The same is true for Eq. ([l.Tp. Therefore we can drop 
all correction terms and we are back in the previous case of D-term lifted supersymmetric 
vacua. 

(2) e<.u<. 
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When we don't impose constraints on the cosmological constant, so we can assume 
\W\ to be generic, we have that the first term in Eq. (4.6) is much smaller than e^^^, while 
cie/u ^ e^/^. The other terms are again much smaller, so ml is dominated by the second 
term. Furthermore {fn^Y' ~ 4|VFp. This is as in the case of pure F-term susy breaking. 

The case with cosmological constant constrained near zero is similar. Since /v? ^ e, 
the l-Fp-part of the potential will dominate over the D-term part, and thus to get near 
zero vacuum energy we should take \W\ F ^ e/u 3> u. This implies that the first term 
in Eq. (|4.6| ) is of order e and can be dropped together with the D-term, bringing us again 
to the situation of pure F-breaking. For similar considerations hold. 

Hence for either case the measure 5'^"^{dV) \ det V"\ becomes 

S^-iF + eN-^d) n (Ai - \Wf ) . ^j;^^ (4.8) 

^ +eiv a) ^1^1^ (4.yj 

with approximated as in the pure F-breaking case. Note furthermore that the delta- 
function forces F to lie approximately in the direction of the eigenvector "^f as u = 
Ai — 2\W\ is small. Finally, since F ~ e/u, we can write 

dXi du dF 
u u F 

and change variables from Ai to F (mapping the domain 2\W\ -|- e ^ Ai <C 2\W\ + e^/^ to 
gi/2 <^ F <^1). This completely reproduces the measure Eq. ( p.4| ) of pure F-breaking vacua 
in this regime. Thus, pure F-breaking vacua with susy breaking scale Mgusy = F ^ e^/^ are 
just slightly perturbed by adding a D-term potential eVo, and their distributions remain 
essentially identical. 

Note that because of the suppression by at least F^,, these vacua can be expected to 
be significantly less numerous than the lifted supersymmetric ones. 
(3) eV2 <u< ei/3: 

This is the intermediate regime, where D-term and F-term effects are of comparable 
size. When e is small, this corresponds to only a small fraction of susy breaking vacua, 
compared to the other two regimes. 

4.3 Summary 

Adding a D-term which becomes small in a region TZ of moduli space will remove all vacua 
located outside that region from the ensemble of vacua with "small" susy breaking scale 
{Msusy n^p), simply because the D-term either destabilizes the vacuum, or it causes 
the susy breaking scale to be too large. Within TZ, it will add new susy breaking vacua 
obtained by lifting susy vacua. 

For the particular case of susy breaking by an anti D3-brane near a conifold point, the 
number of vacua with susy breaking scale Mgusy < will approximately be given by 

M^{Msusy < M,) ~ ^^^^^ (4.10) 



-25- 



where Ngusy (T^) is the number of tachyon-free susy vacua in TZ. Restricting the cosmological 
constant to a small interval of width AA around zero just multiplies this number by AA. 
Apart from these lifted susy vacua, descendants of pure F-term breaking vacua are also 
present and are only slightly perturbed as long as |-Fp is bigger than the added D-term 
potential. Their number can be expected to be roughly 

NpiMsusy < M*) ~ MlNsusy{Tl) (4.11) 

without constraint on A, and 

M'FiMsusy < M,) ~ M}^ AAMsusyin) (4.12) 

if A is constrained near zero. 

Finally, although detailed considerations in the observable sector are beyond our 
scope here, one effect which must be mentioned is the universal — 2|H^p contribution to 
the bosonic mass matrix Eq. ( |2.4D , which will destabilize vacua with light fermions and 
\D\, \W\ » \F\, thus removing a large class of potentially realistic vacua. 



5. Conclusions 



We derived formulae for the distribution of nonsupersymmetric vacua in a general ensemble 
of effective supergravity theories, such as Eq. ( |3.4| ). This distribution is rather similar to 
that for supersymmetric vacua, with certain "correction factors" which we explained. 

We began by reformulating the problem of finding nonsupersymmetric F breaking 
vacua as that of finding eigenvectors of the matrix D^W of second derivatives of the 
superpotential. This makes several features of the problem manifest, most importantly that 
metastable nonsupersymmetric vacua are generic, of number comparable to the number of 
supersymmetric vacua. 



We then argued that the suggestion of |11, ^] for a large power law growth in the 



number of F breaking vacua in models with many moduli is not expected in general. While 
the heuristic argument leading to this suggestion seems sensible, namely that independent 
supersymmetry breakings in independent hidden sectors combine additively and favor high 
scale breaking, our detailed analysis shows that the different F terms will be independently 
distributed only when the different hidden sectors are totally decoupled, which seems un- 
likely to us based on our studies so far. It should be said that this does not address the 
similar argument in [^] in terms of susy breaking by multiple antibranes, and the related 
idea in |22] that multiple independent D terms could lead to the same effect, as the details 
there are quite different, and this possibility remains open. 

We went on to compute the density of nonsupersymmetric vacua, Eq. ( |3.1D , in general 
and in simple examples. Much of the structure of the result comes from the factor det V" 
which is included to normalize the delta function 5{V') to one for each vacuum. Physically, 
this factor is the product of masses squared for all bosons, and this makes precise the 
general expectation that bosons much lighter than the natural scale of the potential (for 
flux potentials, the string scale) are disfavored. 
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Perhaps the simplest summary of the results for pure F breaking is Eq. ( p. 91) , which 
states that the number N of pure F breaking vacua with supersymmetry breaking scale 
Msusy < Af* << Mpi, and with A ~ 0, goes as ~ M^^. This somewhat surprising claim 
is made up of a number of factors, explained in detail in section 3: 

• The most naive expectation would have been a uniform distribution d?F for the 
complex parameter F, leading to AA ~ M^. However, as discussed in section 2.2, the 
equations V = determine the phase of F, leading instead to the generic distribution 
dF and the counting AA ~ M^, for vacua in which the c.c. is not tuned, or in which 
it is tuned by other effects (say by D terms). 

• If we tune the c.c. to zero using |Fp = 3|l^p, we need to take the W distribution 
into account. is a complex variable, and is proportional to the lowest eigenvalue A 
of the matrix Z = D^W (in the sense of Eq. ( p. 17 )). Since Z is complex symmetric, 



this is distributed as d(A^), leading to an extra factor of A ~ ~ F ~ M^. 

The measure factor det V" = Y\ - mf weighs the density of vacua with the product 
of the masses squared of every boson. As discussed in section 2.3, the two bosonic 
partners to the goldstino generically have masses determined by W and F, and this 
leads to a factor ~ Mf . 



Finally, metastability requires a tuning which leads to an additional factor F'^ ~ M, 



* 1 



as discussed in section 2.3. 
In sum, enforcing A ~ has a significant effect on the distribution, apparently favoring high 



scale breaking enough to outweigh factors such as Eq. (L2) and the suggestions discussed 
in section 3.2.2. Note that the arguments only really require A < which is good 

as any more precise requirement would be spoiled by subsequent corrections. 

We can compare this to a rough picture of the "D breaking" vacua. Since the D 
parameters are real, and this type of supersymmetry breaking need not come with light 
scalars, one does not get the measure factors we just discussed, leading to distributions such 
as Eq. (|4.4| ) which would seem to favor low scale breaking. On the other hand, multiple D 
parameters might still lead to power law growth in Msusy, and a full analysis might suggest 
other measure factors, so this remains unresolved. We also do not know whether F or D 
breaking is more common, since we have as yet no explicit ensemble which includes both 
types of vacuum. 

All of these considerations are of course within a subsector of the theory and ignore 
many further corrections. Nevertheless the factors we just listed which entered in our final 
results appear fairly generic, and thus we believe they will be important ingredients in all 
computations of this type. Numerically, they appear as important as the specific physical 
effects focused on in previous work. 

Much more detailed work in examples will be required to judge to what extent the full 
picture is controlled by such generic features. One can take various attitudes about them 
- perhaps they are background, perhaps they are the signal. But they appear to us to be 
very basic properties of the distribution of string/M theory vacua. 
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Of course, we are still only talking about a small part of the full problem of counting 
realistic string/M theory vacua. One might expect equally important selection factors to 
arise at subsequent stages, and it is probably premature to read too much phenomenological 
significance into our results so far. But we would like to suggest that analyses at this level 
of detail, at least for the best understood classes of vacua (Calabi-Yau compactification of 
lib, F and heterotic theories), for every part of the full problem, might be feasible over the 
next few years, and that this would enable us to make similar statements with some claim 
to significance. 
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